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Precursor effects of the superconducting state caused by d-wave phase-fluctuations 

above T c 
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One of the hallmarks of high-temperature superconductors is a pseudogap regime appearing in 
the underdoped cuprates above the superconducting transition temperature T c . The pseudogap 
continously develops out of the superconducting gap. In addition, high-frequency conductivity ex- 
periments show a superconducting scaling of the optical response in the pseudogap regime, pointing 
towards a superconducting origin of the pseudogap. The phase-fluctuation vortex scenario is fur- 
ther supported by the measurement of an unusually large Nernst signal above T c and the recently 
observed field-enhanced diamagnetism which scales with the Nernst signal. In this paper, we use a 
simple phenomenological model to calculate the paraconductivity and magnetic response caused by 
phase fluctuations of the superconducting order parameter above T c . Our results are in agreement 
with experiments such as the superconducting scaling of the optical response and the spin (or Pauli) 
susceptibility, and further strengthen the idea of a phase-fluctuation origin of the pseudogap. 

PACS numbers: 71.10.Fd, 71.27.+a, 74.72.-h 



I. INTRODUCTION 

A variety of different experiments indicate a suppres- 
sion of low-frequency spectral weight in the underdoped 
cuprates below a characteristic temperature T*, that is 
higher than the superconducting (SC) transition temper- 
ature Tl 1 ! 2 : 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 . The size of this pseudogap tem- 
perature regime scales with the superconducting gap. In 
addition, high-frequency conductivity experiments in un- 
derdoped Bi2Sr2CaCu208+5 (Bi2212)ii have indicated 
a SC scaling behavior of the optical conductivity. Fur- 
thermore, a strongly enhanced Nernst signal was mea- 
sured above T c in underdoped samples of La2- a; Sr 3 ;Cu04 
(LSCO), which is usually associated with the presence of 
vortices in the superconducting state and, therefore, im- 
plies that T c corresponds to a loss of phase rigidity rather 
than a vanishing of the pairing amplitude 1 -. Recently, a 
field-enhanced diamagnetism was found in the pseudo- 
gap state of Bi2212, which scales exactly like the Nernst 
signal 12 and further supports the vortex description for 
the loss of phase coherence at T c . 

All these experiments point towards a picture, where 
the pseudogap arises from phase fluctuations of the su- 
perconducting gapi 3 -. In previous workia, we provided 
a detailed numerical solution of a minimal model which, 
however, contains key ideas of the phase-fluctuation sce- 
nario: In this scenario, below a "mean-field" tempera- 
ture Tmf = T* , a c? a ,2_ l) 2-wave gap amplitude is as- 
sumed to develop. However, the SC transition is sup- 
pressed to a considerably lower temperature T c by phase 
fluctuations 13 . In the intermediate temperature regime 
between T* and T c , the phase fluctuations of the gap 
give rise to the pseudogap phenomena. In this previ- 
ous work, it was shown that a two-dimensional BCS-like 
Hamiltonian with a d x 2_ y 2-w&ve gap and phase fluctua- 
tions, which were treated by a Monte-Carlo (MC) simu- 
lation of an XY model, yields results which compare very 



well with scanning tunneling measurements over a wide 
temperature rang o 9 : 14 . Furthermore, this phenomeno- 
logical phase fluctuation model was also able to explain 
the "violation" of the in-plane optical integral in under- 
doped Bi 2 Sr2CaCu20 8 +<5 (Bi2212)i£ and the character- 
istic change of the quasiparticle dispersion observed in 
crossing T c in underdoped Bi2212 samples^. 

Here, we consider in detail two other experiments, 
which critically probe the phase-fluctuation scenario. 
The first one is the optical response. Our calculations 
were, in particular, motivated by optical studies of the 
Berkeley group^i. In this seminal work high-frequency 
optical response, i. e. the complex frequency dependent 
conductivity, was used to capture the short-time scale dy- 
namics of the phase correlations. Indeed, the pseudogap 
conjecture had early on led to an intensive search for nor- 
mal state "remnants" of the SC state, such as the infinite 
d. c. conductivity. However, it was not before this Corson 
et. al. work that it was realized that the short phase- 
correlation times require a high-frequency optical probe 
to make the conductivity enhancement over normal-state 
values observable. 

We give a detailed analysis of the temperature and 
frequency dependence of the optical conductivity <j(lu) — 
<Ji{w) + i(J2(uj) within our phenomenological phase fluc- 
tuation model, which, for completeness, is reintroduced 
in section II. In section IIIII we compare the calculated 
paraconductivity due to cZ-wave phase fluctuations above 
T c with the high-frequency conductivity experiments of 
Ref. ^3- Salient features of these experiments are re- 
produced in (t(oj) above T c : at sufficiently high frequen- 
cies, the phase-fluctuating state becomes indistinguish- 
able from the SC state. In this latter case 02 (w), i. e. the 
imaginary part, shows the usual BCS-like behavior falling 
off as l/u>, with a prefactor being proportional to the 
phase-stiffness energy. This generates in our calculations 
a scaling behavior of <72(w,T), in agreement with exper- 
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iment, in that all ujg2{oj,T) curves collapse onto a sin- 
gle curve at the Kosterlitz-Thouless (KT)-temperature 
T c = Tkt- 

Since, in the phase fluctuation scenario, the pseudogap 
is due to phase fluctuations of the SC order parameter, 
one might ask, how much of the characteristic magnetic 
properties of a superconductor are still observable in the 
pseudogap state of the underdoped cuprates. 

These magnetic precursor effects of the ideal supercon- 
ducting state are twofold. Firstly, one would expect a 
form of fluctuating diamagnctism which partially screens 
an applied magnetic field already above T c . Secondly, one 
would also expect that the paramagnetic spin (or Pauli) 
susceptibility is reduced in a characteristic way below T* 
due to the formation of incoherent singlet-pairs, accom- 
panied by the reduction of spectral weight at the Fermi 
surface. 

In experiment, the fluctuation-induced diamagnetism 
is only important close to T c . Since, at higher tem- 
peratures, the paramagnetic spin response is typically 
one-order-of-magnitude larger than the diamagnetic re- 
sponse, we concentrate here on the spin susceptibility Xs- 
This latter quantity dominates the T-dependence of the 
paramagnetic susceptibility in the experiments. 

In section I1VI it is shown that the Pauli spin sus- 
ceptibility extracted from our phase-fluctuation model 
can qualitatively account for the observed paramagnetic 
properties over a wide range of temperatures in the pseu- 
dogap regime. Section V, finally, presents a summary of 
our results. 



II. PHASE-FLUCTUATION MODEL 

We consider the Hamiltonianii 

H = H + H U (1) 

where Hq is the usual tight-binding operator of non- 
interacting electrons on a two-dimensional (2D) square 
lattice, i. e. 

H ° = - t E ( c L c ^+4. c - ) - /i E^- ( 2 ) 

,CT 1,(7 

Here, ct (c? CT ) creates (annihilates) an electron of spin a 

on the i th site of the (2D) square lattice and n? = ct c? 
is the number operator, t denotes an effective nearest- 
neighbor hopping-term and u, is the chemical potential. 
The (ij) sum is over nearest-neighbor sites of the 2D 
square lattice. For simplicity of our minimal model, we 
set longer-ranged hoppings t' equal to zero. All calcula- 
tions were done for 10 % doping, which corresponds to a 
typical underdoped cuprate situation. 

The second part of the Hamiltonian Hi contains a 
BCS-like <i-wave interaction, which is given by 

tfi = -sE( A ^< A ^>+ A M A ^)< ( 3 ) 

Is 



with 5 connecting site i to its nearest-neighbor sites. The 
coupling constant g stands for the strength of the effec- 
tive next-neighbor d x 2_ y 2-w&ve pairing-interaction. The 
origin of this pairing interaction is not essential for the 
further calculation. It can be either of purely electronic 
origin, like spin fluctuations, or phonon mediated. The 
only important point is, that there exists an effective 
pairing interaction, that produces a finite d x 2_ y 2-w&ve 
gap amplitude as one goes below a certain temperature 
T*. In contrast to conventional BCS theory, we consider 
the pairing-field amplitude not as a constant real number, 
but rather as a complex quantity 

(At-) = -L(ct ct - - etc! - ) = Ae«*«, (4) 

with a constant magnitude A and a Euctuating bond- 
phase field In order to get a description, where the 
center of mass phases of the Cooper pairs are the only rel- 
evant degrees of freedom^, the d x 2_ J( 2-wave bond-phase 
field is written in the following way 

$ _ I (<p*+ <Pj + g)/2 for Jinx-direction 

1 (ip^ + ip* + g)/2 + 7r for 6 in y-direction, 

where ip? is the center of mass phase of a Cooper pair 
localized at lattice site i. 

To account for the proximity to the Mott insulating 
state and thus the low superfluid density, we perform a 
quenched average over all possible phase configurations 
with the statistical weight given by the classical XY free 
energy 

F[pi] = -J^2cos((fi - (fj). (6) 

(«) 

Here, the phase stiffness J determines the Kosterlitz- 
Thouless transition temperature Tkt to a quasi phase- 
ordered state, a temperature, which we take as 
Finally, we set T c ss \T* , where we have the STM exper- 
iments of Ref. in mind, which corresponds to a large 
pseudogap regime . 

In Ref. [ij, the justification for using a minimal model 
and, in particular, the use of a classical XY interaction in 
Eq. ©, has been carefully discussed. Our belief, that this 
model contains key ideas of the cuprate phase-fluctuation 
scenario, has been substantiated there by a detailed nu- 
merical solution and the fact that characteristic and cru- 
cial features of STM studies^ could be reproduced over a 
wide temperature range. 

Recent work on a BCS-Hamiltonian with classical 
phase fluctuations can also be found in Ref. 

III. PARACONDUCTIVITY 

In a conventional BCS superconductor, superconduct- 
ing (SC) fluctuations with short-range phase coherence 
typically survive no more than IK above T c . Within a 



3 



phase fluctuation scenario for the underdoped cuprates, 
one would expect that pairing remains over a wide tem- 
perature range above T c , together with phase correlations 
which arc of finite range in space and time^. Hence, al- 
though the system is in the normal state, signatures of 
the ideal SC state should still be observable considerably 
above T c , if the experiments probe short enough time 
scales. A likely candidate to observe these effects are 
high-frequency conductivity experiments. 

Indeed, microwave conductivity experiments on under- 
doped Bi2213il were able to track the phase-correlation 
time in the normal state up to 25K above T c . These 
experiments show, that the SC transition is "smeared 
out" over a considerable temperature range, when viewed 
at high-enough frequencies (short-enough time scales). 
The imaginary part of the conductivity finally disappears 
more than 25K above T c ~ 74K and shows a supercon- 
ducting scaling behavior already above T c . The real part 
of the conductivity displays a characteristic peak near T c 
at finite frequencies, on top of a "background conduc- 
tivity" of normal conducting electrons. This peak was 
interpreted by the authors of Ref. ^] as signature of the 
partially phase coherent electrons. 

All experiments were carried out with frequencies of 
a few hundred GHz, which corresponds to ui < O.Olt. 
In order to compare these experiments with our phase 
fluctuation model, some compromise has to be made, 
since finite-size effects become important below u>i ow — 
O.QSt. Therefore, we have chosen a set of frequencies 
u> = {O.lt, 0.2i, 0.3i, OAt}, which are, on the one hand, 
much larger than u>i ow , and on the other hand much 
smaller than the relevant energy scale, i. e. the SC gap 
size in our model (A sc = l.Ot). Thus, in the following, 
one always should keep in mind that the frequency closest 
to experiment is u) = O.lt. 

The optical conductivity from our phase-fluctuation 
model, i. e.— 

ft \ 2 (M -T> c (k,u + ir)) 
a xx (k,u)=e — — , (7) 

l(UJ + 11]) 

was obtained by averaging the Matsubara current- 
current correlation function 

V c (lr, f r<) = -Tr{p G T r \£(T, t) j£(f, r')]}, (8) 

and the operator for the local kinetic energy in in- 
direction 

k x (l) = -t V(ct c r +ct cr,- ) (9) 

a 

over all possible phase configurations using Eq. JJjJ. Here, 
Pg is the usual statistical operator, and 

fjT) =itV(ct c T -cl cr- ) (10) 

cr 

denotes the paramagnetic current density operator in in- 
direction. 
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FIG. 1: Imaginary part of the low- frequency optical conduc- 
tivity 0-2(0), T), as a function of temperature for different fre- 
quencies ui < A 3C (A ac = l.Ot). The dotted vertical line 
indicates T c = Tkt- Note, that the change of <72(w,T) at 
the superconducting transition T c = Tkt becomes less pro- 
nounced at higher frequencies and is smeared out over a finite 
range of temperatures T > T c = Tkt- 







1 1 1 1 1 

jnJC^"*—* A ,v 






^3-/""* e .. 






G-O(0=0.1t 
Q-B (O=0.2t 

O-&(O=0.3t 

A^(B=0.4t 

i 


I.I. 



1.0 2.0 3.0 4.0 

T/Tkt 

FIG. 2: Scaling behavior of the imaginary part of the op- 
tical conductivity U2(lo,T). In the superconducting state 
<72(w,T) ~ thus all uj (T2(w,T) curves should collapse 
onto a single curve. Note, that the higher the frequency, the 
earlier this collapse starts in the normal state for tempera- 
tures T > T c = Tkt- The logarithmic scale was chosen for 
a better comparison with the experimental data of Ref. ITU 
reproduced in Fig. 3, below. 



Fig.^plots our results for the imaginary part of the op- 
tical conductivity,!. e. 02(0;, T), as a function of tempera- 
ture for different frequencies u> < A sc . Exactly as in the 
experiments of Ref. Illl the superconductivity induced 
change in the imaginary part extends over a large temper- 
ature interval T c < T < 2 T c into the normal state, with 
the strongest increase below T ~ 1.5 T c . The change of 
<72(w,T) is more pronounced at lower frequencies. How- 
ever, in our model the imaginary part does not go to zero 
in the normal state but rather converges towards a finite 
value. This is due to the missing quasiparticle scatter- 
ing in the normal state in addition to phase-fluctuation 
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FIG. 3: Experimental results for the dynamic (frequency de- 
pendent) phase- stiffness temperature Tg ~ ui azjij.T) in the 
pseudogap state of underdoped Bi2212 from Ref. HH Data 
are shown for two samples, one with T c — 33K (left side) and 
the other with T c — 74K (right side). The three curves for 
each sample correspond to measurement frequencies of 100, 
200 and 600 GHz. (For details, see text.) 



caused scattering. In a "real" metal, above T c , the finite 
electronic scattering rate, strongly suppresses the imag- 
inary part (72 (w) at low frequencies, so that it can be 
neglected. This effect can phenomenologically be taken 
into account by replacing the infinitesimal damping fac- 
tor r\ in Eq. Q by a finite marginal-Fermi-liquid (MFL) 
damping factor 19,20 . Since the inclusion of such a MFL 
scattering-rate further reduces <72(w,T) above T ~ 2T C , 
however without changing its properties at the SC transi- 
tion (T < 2 T c ) qualitatively, we refrain from burdening 
our simple model with it. 

In Fig. |21 the re-scaled imaginary part of the optical 
conductivity w a%{u), T) is displayed as a function of tem- 
perature for the same frequencies to < A sc , as before. In 
the superconducting state, all re-scaled curves should col- 
lapse onto a single curve, as has been beautifully demon- 
strated in Ref.[Tl| (Fig. 2 of that work, which corresponds 
to our Figgis reproduced for comparison in Fig- EJ ■ One 
can clearly see from our results that this collapse already 
begins in the normal state above T c = Tkt, starting 
with the highest frequencies, exactly as in Ref. 0. Also 
here, the inclusion of an additional MFL scattering-rate 
reduces the optical conductivity curves above T ~ 2T C , 
but it does not change the SC high-frequency scaling- 
behavior below T < 1.5 T c significantly^. 



IV. MAGNETIC PROPERTIES 

In this section, we study the precursor effects of the 
SC state caused by d-wave phase fluctuations above T c 
appearing in the paramagnetic response. The spin (or 
Pauli) susceptibility can be obtained from the coupling 
of the magnetic field B to the spin degrees of freedom. 
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FIG. 4: Spin susceptibility Xspm in units of g 2 fi%/t (h — c = 
a = 1) as a function of temperature for a constant d x i_ y 2- 
wave gap with phase fluctuations ((n) = 0.9). For compar- 
ison, we also show the temperature dependence of the spin 
susceptibility Xspin for a BCS-superconductor with a constant 
gap and for the tight-binding model (A = 0). 



With B || z, one getsSi 

Xs P in(q,u) = - 



N 2 Ml 1 n ^ x 



(11) 



where D zz (q,uj) is the Fourier transform of the spin-spin 
correlation function in z-direction, i. e. 

D zz (x,t, x',t') = -iTr{p G [s z {x,t),s z (x',t')}}0(t-t'). 

The static uniform susceptibility is defined by Xspin = 
Xspin(q — * 0, lu = 0). For a non-interacting system with a 
general dispersion e(k), the static uniform spin (or Pauli) 
susceptibility is given by 



Xs P in(q -> 0,w = 0) 



1 

V 



2 2 
9 (J-B 



E 

k 



df(e(k) - /x) 
de(k) 



(13) 

with f(e(k) — /i) being the Fermi function. For T — > 0, 
the static uniform susceptibility is proportional to the 
density of states at the Fermi surface. 

Fig. 2] displays the spin susceptibility \spin as a func- 
tion of temperature for a constant d x 2_ y 2-w&ve gap with 
phase fluctuations at finite doping ((n) = 0.9). In Sec. [I] 
we have mentioned already that the fluctuation-induced 
diamagnetism is only important below T ~ 1.5 T c and 
that, at higher temperatures, XDia is at least two or- 
ders of magnitude smaller than the paramagnetic current 
response of the tight-binding electrons. Moreover, the 
experimentally observed diamagnetic susceptibility close 
to T c can be perfectly fitted by an exponential function, 
which also indicates that the fluctuating diamagnetism is 
exponentially suppressed at higher temperatures^. Fur- 
thermore, one obtains that the paramagnetic spin re- 
sponse is about 5 times larger than the paramagnetic 
current response of the tight-binding electrons. In addi- 
tion, the magnetic current response is only weakly tern- 
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FIG. 5: Spin susceptibility Xspin in units of g 2 n 2 B /t (h = 
c — a = 1) as a function of temperature for a BCS- 
temperature-dependent d a ,2_ y 2-wave gap with phase fluctu- 
ations ((n) = 0.9). For comparison, we also show the temper- 
ature dependence of the spin susceptibility Xspin f° r a BCS- 
superconductor and for the tight-binding model (A = 0). 
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FIG. 6: Experimental results for the c-axis NMR Knight-shift 
in underdoped, optimally doped and overdoped Bi2212 from 
Ref. I24L The temperature variation of the NMR Knight-shift 
scales linearly with the macroscopic magnetic susceptibility 
Xspin- The anomalies indicated by arrows are discussed in 
the text. 



perature dependent. Hence, we expect that the temper- 
ature dependence of the experimentally observed param- 
agnetic susceptibility is dominated by the spin suscepti- 
bility Xspin over almost the the entire pseudogap phased 
In Fig. E| 

Xspin is plotted for our phase-fluctuation 
model and, additionally, for the nearest-neighbor tight- 
binding model and for a BCS-model with a constant 
superconducting gap. The phase-fluctuation model ex- 
hibits a characteristic temperature dependence of the 
spin susceptibility, which differs qualitatively from the 
tight-binding model as well as from the BCS-model. 
This becomes even more clear for a BCS-temperature- 
dependent gap, as shown in Fig. |SJ From Fig. [S] we can 



also infer, that a temperature dependent pairing-gap only 
slightly modifies the temperature dependence of Xspin m 
our phase-fluctuation model. Xspin always slightly de- 
creases (nearly linear) below T* = T C MF and, then, dis- 
plays a characteristic downward bending at T ~ 2Tkt- 
These results are very similar to the experimentally ob- 
served temperature dependence of the NMR Knight-shift 
and the magnetic susceptibility in the pseudogap state 
of various underdoped high-T c compounds^2SiSS42i. For 
comparison, we reproduce in Fig. the T-dependence of 
the Knight-shift experimental data from Ref.l24lin under- 
doped and overdoped Bi2212. In the pseudogap state, the 
temperature variation of the NMR Knight-shift scales lin- 
early with the macroscopic magnetic susceptibilitj* 2 ^. In 
underdoped Bi 2 Sr 2 CaCu20g + 5 (Bi2212) single-crystals, 
two characteristic temperature scales can be identified in 
the temperature dependence of the Knight-shift 24 (sec 
Fig. EJ. A higher temperature T m K, where the Knight- 
shift starts to decrease from the nearly constant high- 
temperature value, and a lower temperature T% > T c , 
where it starts to decrease very steeplj* 2 ^. These two tem- 
peratures are, however, difficult to define exactly. This 
can be seen from the scaled Knight-shift date of YBCO 
shown in Ref. |2{| and also from our numerical results dis- 
played in Fig.0]and Fig. [SI since the change as a function 
of temperature is continuous. 



V. SUMMARY AND CONCLUSION 

In the high-T c superconductors, a pseudogap opens at 
at temperature T* that is higher than the SC transi- 
tion temperature T c . There is a rather general consensus 
that the pseudogap state is crucial for understanding the 
microscopic pairing mechanism, however key issues are 
unsettled. One of these is whether it is a precursor state 
for SC, sharing common features such as pair formation 
or whether it is even antagonistic to SC. 

At present, there is not yet agreement as to which of 
these possibilities is correct. In part, this is because of the 
difficulty in determining the experimental consequences 
of the various theoretical proposals. This is exactly what 
we seek to improve here. 

We have provided a numerical solution of a simpli- 
fied model which, nevertheless, contains key ideas of the 
phase-fluctuation pseudogap scenario. Here, the center- 
of-mass pair-phase fluctuations of a BCS d-wave model 
were determined from a classical 2D X Y-action by means 
of a Monte-Carlo simulation. Earlier work concerned 
with the one-particle excitations have already been found 
to reproduce salient features of recent STM studies of 
Bi2212 and Bi2201— and to explain the change of the 
quasiparticle dispersion in crossing T c in underdoped 
cupratesi^. Here, we concentrate on experiments related 
to two-particle correlation functions, i. e. the optical and 
magnetic response functions. They are sensitive, diag- 
nostic tools to search for normal-state remnants of the in- 
finite d. c. conductivity and perfect diamagnetism above 



6 



In the first part of the present work, we have studied 
the frequency-dependent complex conductivity. Our mo- 
tivation there stems partly from the work of the Berke- 
ley groupii, which very convincingly demonstrated that 
the consequences of partial phase coherence can only be 
resolved by high-frequency optical experiments, probing 
the short-time dynamics of the phase correlations. 

A detailed analysis of the temperature and frequency 
dependence of the optical conductivity cr(u>) = <Ji(u>) + 
ia2(to) revealed a superconducting scaling of <t 2 (o;), 
which starts already above T c , exactly as observed in 
high-frequency microwave conductivity experiments on 
Bi2213±i. 

Secondly, we calculated the magnetic susceptibility. 
The temperature dependence of the uniform static 
magnetic susceptibility is dominated by the Pauli spin 



susceptibility, which displayed a very characteristic 
temperature dependence, independent of the details 
of the gap function used in our model. This temper- 
ature dependence is qualitatively very similar to the 
experimentally observed change of the Knight-shift as a 
function of temperature in undcrdoped Bi2212. 
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